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^ Abstract 

^ The question of adaptive mesh generation for approximation by sphnes has been studied for a 

^ number of years by various authors. The results have numerous apphcations in computational and 
y—i discrete geometry, computer aided geometric design, finite element methods for numerical solutions 
^ of partial differential equations, image processing, and mesh generation for computer graphics, 

among others. In this paper we will investigate the questions of adaptive approximation of 
^ functions with arbitrary but fixed throughout the domain signature by multilinear splines. In 
fH particular, we will study the asymptotic behavior of the optimal error of the weighted uniform 
approximation by interpolating and quasi interpolating multilinear splines. 
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2^ 1 Introduction. 

Let the domain D be the unit cube [0, l]'^ C R'^, d eN and d>2. However, any bounded connected 
^ region that can be represented as a finite union of cubes can be treated analogously. 
2 Let Loo(-D) be the standard space of essentially bounded measurable functions defined on D with 

^ the usual sup-norm || ■ ||oo- Given a positive continuous function Q{x.) on D, define the weighted 
^ uniform norm || • ||oo,n of / G Loo(-D) as 

>< \\f\\^,n := sup{ |/(x)|Q(x) : x e L> } = \\fn\\^. 

Let Div = {Ri}fL^ be an arbitrary partition of the domain D into parallelepipeds (or boxes) 
Ri with sides parallel to coordinate axes. We shall refer to such a partition as a box partition. 

Let P* be the space of polynomials linear in each of its d variables. Given a function / e L^{D) 
and a partition Djv we shall consider the multilinear spline Dat) defined as follows: 

1. On the interior of each box i = 1, . . . , A?", we define D^v) to be the polynomial from 

which interpolates / at the vertices of 

2. For every point x in the union of the boundaries of boxes Ri, i — 1, . . . , N we define the value 
of Div) to be the average value at x of all polynomial interpolants on the boxes whose 
boundaries contain x. 
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Observe that, by construction, splines s{f, O^) interpolate the function / at the vertices of the 
partition Dtv- However, in general, they are not necessarily continuous. 

We shall call a sequence of partitions {DAr}^^]^ admissible if it satisfies the condition 

sup A^^ max diam(i?) < oo. (1) 

TV -Renjv 

Throughout this paper we shall consider only admissible box partitions. 
The quantity 

R^(/):=inf||/-s(/,n^)|U,n (2) 

Ujv 

where inf is taken over all admissible box partitions Dat of the domain D into boxes we shall 
call the optimal error of interpolation. An explicit form and the exact value of RAr(/), as well as 
the explicit construction of the optimal partition for every particular function /, can be found only 
in trivial situations. It was shown by Below, De Loera, and Richter-Gebert in 2000 [3] that it is 
not possible to construct an adaptive algorithm for optimal mesh (partition) generation that runs 
in polynomial time. 

That is why it is interesting to study the asymptotics of the optimal error Rjv(/) as A^ — oo 
for a given function / G and to construct an asymptotically optimal sequence of box partitions, 
i. e. a sequence of partitions {D^}^^! of D such that 

Note that the problem formulated above is interesting for functions of arbitrary smoothness as 
well as for various classes of splines (for instance, for splines of higher order, interpolating splines, 
best approximating splines, best one-sided approximating splines, etc.). In the univariate case 
general questions of this type have been investigated by many authors. The results obtained in this 
case are more or less complete and have numerous applications (see, for example, [23]). 

Fewer results are known in the multivariate case. The classical statement of L. Fejes Toth 
indicated in [17j on approximation of convex bodies by inscribed polytopes in Hausdorff metric can 
be considered as the first result in this direction. Gruber [18] proved this result and generalized it 
to the arbitrary dimension. We [2] proved similar result for the weighted uniform norm. Related 
interesting results on approximation of convex bodies by various polytopes have been obtained by 
Boroczky [7], Boroczky and Ludwig [S]. 

In [21] Nadler solved the problem of asymptotically optimal choice of a sequence of triangulations 
for approximation of functions by piecewise linear functions of best L2-approximation. 

D'Azevedo and Simpson [Hj studied the question of triangulating a given set of vertices for 
interpolation of a convex quadratic surface by piecewise linear functions. They showed that the 
Delaunay triangulation will be optimal for the error in the L^o norm. For the error in Lp norm 
this fact was proved by Rippa [26] . Chen, Sun, and Xu [TT] generalized this result to arbitrary 
dimensions. 

Later D'Azevedo [T3| obtained local error estimates for functions with both positive and negative 
curvature. The same estimates were later obtained by Pottmann, Hamann et al [25] who studied 
the problem of optimally triangulating the plane for approximating quadratic functions by piecewise 
linear functions and suggested some algorithms for constructing function dependent triangulations 
of the whole domain. 
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Huang [21] , and Huang, Xu, and Sun considered the problem of variational mesh adaptation 
in the numerical solutions of partial differential equations and obtained asymptotic bounds on the 
interpolation error estimates in L2 for adaptive meshes that satisfy regularity and equidistribution 
conditions. 

All above mentioned results are for the case of approximation by linear splines. Natural domain 
partitions in this case are simplices. However, in applications where preferred directions exist, box 
partitions (or generalized rectangular partitions) are sometimes more convenient. The only known 
to us result for box partitions is the estimate of the uniform error of interpolation on rectangular 
partition by bivariate splines linear in each variable which is due to D'Azevedo who obtained 
the local error estimates in this case. 

Therefore, we think it is an interesting problem to obtain results similar to the above mentioned 
for the class of splines linear in each variable (multilinear splines) defined over box partitions. 

In this paper we shall construct an asymptotically optimal sequence {D^}^^^ of box partitions 
and determine the exact asymptotic value of RAr(/) for a function / G C'^{D) with fixed signature 
(with fixed number of positive and negative second derivatives) throughout the domain. Even 
though signature is fixed, the results presented are richer than what has been proved for interpolating 
by linear splines on M'^, d > 2. In the latter case, while some results exist for positive definite 
functions, no results exist in the case of functions of another signature. 

The approach we shall undertake is similar to the one used for studying the asymptotics of the 
error of approximation by linear splines and is as follows: we first take an intermediate approxima- 
tion of the given function / by a quadratic polynomial and then find the error of approximating the 
quadratic by multilinear splines. This last problem is solved for quadratic functions of arbitrary 
fixed signature. The result is then used to give the error of approximating a function with (fixed 
throughout the domain) arbitrary signature. Even though it has not being done in the current text, 
it is rather clear how to proceed and extend the obtained results to approximate an arbitrary 
function. 

In addition, we show that the corresponding multilinear splines {■5(/, □^)}5^=i can be con- 
structed so that they will be discontinuous only along small number (in comparison with the total 
number) of lines. 

If we do not require interpolation at every vertex of a partition we can construct an asymp- 
totically optimal sequence of admissible partitions {□Ar}^=i and a sequence of continuous splines 
{s(/, 0^)}^^^ which interpolate / at all but o{N) vertices of the partition as — )■ 00. We shall 
refer to such splines as quasi interpolating splines. 



2 Main results and ideas of proofs. 



For / G C2 



{D) denote 




d 



(4) 



i=l 



Let us consider for any < A; < c? the following class of functions 



cm := { 



/ G C^{D) : Vx G 



(x) > 0, l<i<k, and 



(x) < 0, k <i <d 



dxf 



dxj 
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Sometimes, we shall say that functions from Cl{D) have signature {k,d). In the case when k = 
OT k = d we shall say that functions are positive definite. 
In addition, for fc, G N U {0} set 



lk's{d-ky-'^, 0<k<d 

k = d OT k = 0. 



The next theorem contains the main result of this paper. 
Theorem 1 For any < k < d and f G Cl{D) 



, 2 . 
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hm iV3R^(/)=7Mf / |//(/;x)|5l](x)5dxy. (6) 

Furthermore, there exists a sequence of admissible box partitions {\^n}^=i and a sequence of 
continuous quasi interpolating splines Dat) such that 

^.^ ||/-g(/,D^)||oo,n ^ ^ 
N^co R^(/) 

Remark. The sequence of splines constructed in the proof of the upper bound in ^ will possess 
the following nice property: for each N the constructed spline will be discontinuous along only a 
small number, i.e. o{N) of faces as — )■ oo (see Section 5.2). 

Let us describe the idea of the proof of the estimate from above. It consists of finding an 
appropriate sequence of "good" partitions of D. This is done in the following way: 

1. Divide D into a number (which is small in comparison with A^) of equal subregions . 
On each Df , instead of /, consider the quadratic part of its Taylor polynomial taken at the 
center hj of Df (call it Pl{f,^)). The error of this intermediate approximation is given in 
Lemma [H 

The choice of m% is governed by the following two reasons: there should be few original 
subregions in comparison with A^, but their size should be small enough to provide the small 
enough error of intermediate approximation of / by P^if,^.). 

2. We find the parameters of the appropriate partition of by minimizing the error of multi- 
linear interpolation of P2{f, x) on Df, i = 1, . . . , mf^. We choose the size of R in such a way 



m 



N 



that the overall number nf of elements of partition used for Df^ is such that the sum 

i=l 

is approximately A^, and the errors of interpolation on each Df are approximately equal. 



3. The final partition of D is the union of partitions of each region Df, i = 1, . . . , m^. 



4. We show that the sequence of partitions which is optimal for the intermediate approximant 
(piecewise quadratic function P2{f, x)) will be asymptotically optimal for the original function 

/■ 
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Having constructed a partition for the fixed A^, we interpolate / at the vertices of this partition. 
This will produce a multilinear spline which will be discontinuous along "small" number of edges 
in the partition. Repeating the construction for every N, we shall obtain a sequence of partitions 
and therefore a sequence of interpolating multilinear splines which will be asymptotically optimal. 

If we "give up" the interpolation at some points ( "small" amount of them) of the sequence of 
partitions for the sake of having a sequence of continuous multilinear splines, then we shall refine 
the obtained on each step partition and then "glue" splines on the neighboring elements together. 



The resulting continuous spline we shall call quasi interpolating spline (see Section 5.2 for detailed 
construction) . 

3 Auxiliary statements. 

The proofs of the following auxiliary statements are straightforward. Similar statements have been 
proved, for instance, in [31 12]. 
Set for f E L^{D) , D C 

a;(/,5):=sup{|/(x)-/(x')|: |x-x'|<5, x,x'gD}, 6>0, (7) 

where Ixl := max IxA for x G M'^. Set for f G C'^(D) 

u*{f,6) := max {uj{f^^^ 6)}, (8) 

l<i,j<d 

where fx^xj denotes mixed derivative of / with respect to variables Xj and Xi, i,j = l,...,d. 

Lemma 1 Let f G C^{D) and -P2(x) denote its quadratic Taylor polynomial at the center xq of a 
cube Dfi d D inW^ with side length equal to h. Then we have the following estimate: 

|/(x) - P2(x)| < ^ u* (f, |V X G D,, (9) 



2 \2 J V 2^ 
where u*{f,t) is defined in 

For a fixed a G M'^ and an arbitrary box R denote 

i? + a={x + a, X G -R}. 
Lemma 2 For the given quadratic function 

d 

g(x) = ^A,x^, (10) 

i=l 

any box R, and a G M*^ errors (in any Lp norm) of multilinear interpolation of Q{x) at the vertices 
of R and R + a. are equal. 



i=l 



Lemma 3 The interpolant of the quadratic function (10) on the d-dimensional box R'^ := ]^[— /ij, hi 
is a constant function 

d 
i=l 

Proofs of last two statements are simple linear algebra exercises. 
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4 Interpolation of quadratic functions with arbitrary sig- 
nature. 

Let the quadratic form 

i=l i=k+l 

for < k < d he given. 



Lemma 4 Lei Q be the quadratic function of form (11), let := ]^[— /ij, /ij], and let P be the 

I i=l 

unique polynomial from interpolating Q at the vertices of R^. Then 

\\Q - Ph^iR") = max{hl + --- + hi hl^, + ... + hl}. (12) 

Proof: We shall proceed by induction on the number of variables. To prove the basis of 
induction we need to consider two cases d = 2, k = 1 and d = 2, k = 2. 
Case 1 . Let the quadratic form 

Q{x,y) =x'^ -y"^ 

and an arbitrary rectangle R := [—hi, hi] x [— /i2, h2] be given. 

By Lemma [sj the bilinear interpolant to the function Q(x,y) on the rectangle R, denoted by 
TQ,R(x,y), is a constant equal to 

TQ^R{x,y) = hl- hi. 

Observe that, because of the symmetry, the error in the uniform norm on R is the same as the error 
on [0,/i2] X [0,/ii]. 

Denote the difference between the function Q{x,y) and the interpolant TQ^ji{x,y) by 

6{x,y) ■.= x^-y'-hl + hl (13) 

Clearly, the point (0,0) is the only critical critical point in [0, /i2] x [0, /ii] of this function. The 
value of the difference (13) at this point is 

mo)\ = \hj-hi\. 

In addition, observe that on the boundary of [0, h2] x [0, hi] we have 

5{x,h2) = — h\, and 5{hi,y) = y'^ — h\ 
and, hence, maximal values of the difference are 

|5(0,/i2)| = /i? and \5{hi,Q)\ = hi 
Therefore, the error in the uniform norm can be rewritten as follows: 

max{|/i^ — hi, hi hi} = max{/i^, hi}, (14) 
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and the statement of the Lemma is proved in the Case 1. 

Case 2 . Similarly to the Case 1 we can show that the error of interpolation of the quadratic 
form 

Q{x,y) = x^ + y'^ 

on an arbitrary rectangle R = [—hi, hi] x [—h2, h2] by multilinear polynomial which in this case is 
going to have a form 

Tq,b{x, y) = hj + hj. 

In addition, note that the error in the uniform norm on the rectangle R is the same as the error on 

[0,/l2] X [0,hi]. 

Denote the difference between function Q{x,y) and interpolant TQ ji{x,y) by 

6{x,y) ■.= x^ + y^-hl-hl (15) 

Clearly, the point (0, 0) is the only critical point of this function inside R. The value of difference 
(15) at this point is 

5(0,0) = hl + hl 

In addition, observe that on the boundary of [0, /i2] x [0, /ii] we have 

5{x,h2) = x^ — h\, and 5{hi,y) = y'^ — h\ 

and, hence, the maximal values are 

|5(0,/i2)| = /i? and \5{hi,Q)\ = hi 

Therefore, in Case 2 the error is equal to 

m&^{h\ + hl,hl,hl} = hl + hl. (16) 

This completes the proof of the basis of induction. 

Next let us consider form (11) with signature (fc, d — k) when < k < d. As before we can see 
that the interpolant to (11) on i? is the constant 

k d 



1=1 i=k+l 

Denote the difference between the function and the interpolant by 

5(x) :=g(x)-TQ,K(x). 
Let us investigate critical points of this function. For brevity of computation, let us denote 

i=l j=ld¥^i i=k+l j=k+l,j^i 

As before, the center is the only critical point of 5(x) inside R. The value of the difference at the 
center is 

5{Q) = \Sk-S,\. (17) 
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Let us consider the error on the boundary. 



On the face Xi = hi in the case when i < k form (11) becomes 

k d 

E ^? - E ^1 (18) 

J=ljy« j=k+l 

and the error by hypothesis of induction is 

mB.^[s^\Sd]. (19) 



Similarly, on the face Xi = hi in the case when i > k form (11) becomes 

E-?- E -I (20) 

j=l j=k+l,j^i 

and the error by hypothesis of induction is 

max{5fc,5«}. (21) 

Therefore, the global error is 

A = max < \Sk — Sd\ ■, Ton&x max < Su\ Sd\ , max max < 5*^, S'l*'' 

l<j<A: { ) k+l<i<d { . 



max < max max < S'[*\ \ , max max < Sh, Sf' \ \ 
[i<i<fc I J fc+i<i<d I ) ] 

max "I max <| max 5"^*^ 5"^ > , max < S'^, max S^f^ 



l<i<k I I k+l<i<d 



= max < max 5'[*'*, S'fc, S'd, max M*^ 
|^i<i<fc fc+i<i<d 

= max {Sk, Sd} ■ 

The lemma is proved. □ 

Next we shall compute the minimal value of the error A for the quadratic form (11) with 
signature {k, d — k). Denote by 

A:=mm|x:^l E ^4 ' ^''^ 

d 

where min is taken over all hi such that the volume 2"^ JJ^ hi is fixed (= V"). In the next two lemmas 

i=l 



we shall provide the value of A in the case of the quadratic form (11 ) with signature (fc, d — k) when 



k ^ 0,d (Lemma 



5) and in the case of positive definite quadratic form Q(x) = ^^x^ (Lemma 



i=l 



6). 



Lemma 5 The minimal L^q- error of interpolation of quadratic form (11) by a polynomial from 
on all d- dimensional boxes of fixed volume V is 



\ k , , -I fc 2 



(23) 



Proof: Due to Lemma 



we may consider only boxes = ]^[— /^i, hi] centered at the origin. 



i=l 



(k d \ 

^^/i^, h^ j given by expression in (12), we 
j=i j=k+i ) 

minimize the /p-norm of this vector (with an arbitrary p) of (p!2|, i.e. the expression 



k_ \ P 



(24) 



Vi=A:+l 



under assumption that the volume of the box is fixed (= V) and take the value of the minimum 
when p = oo. 

Indeed, if for some set M C M'^ we denote 



|a;*||; = inf = y4„. 



then it is not difficult to check that 



= inf llxl 



lim Ap = Aoo- 



A. 



(25) 



Indeed, obviously A^o < Ap, for any p. On the other hand, 

Ar, = \\x*\\l < \\x*\\l < \\x*\\l + Er,. 

P II p\\''p II llt-p II ll'-OO ' P 



Last two estimates combined imply (25). 

Denote the length of the sides of the box by /ij, i = 1, . . . , (i. The assumption of volume being 
fixed is equivalent to 



2^''\{hl = V\ 



(26) 



i=l 



The standard routine of minimizing ( [24|) leads to the only solution of the minimization problem 



h\ = ■ ■ ■ = h\ =: X, h 



fc+i 



hi =■■ y. 



Taking this into consideration together with assumption (26) which now can be rewritten as 



^kyd-k ^ Y22-2d^ 



we can find x and y: 



In the case p = oo we have 



1 fd-kV^^-^K^ 



X 



'''2 , 1 f d-k\ " 

V-^, and y = -A^ 



Therefore, 



d-k 

1 fd-k\ 2d _i 



^ 2 \ k 



Hence, the minimal value A of the error A is 



1 f d-k\ 2d _i 



~ 1, (d-k 



Vi, j > k. 



(28) 
(29) 



(30) 



□ 



Lemma 6 The minimal L^o-e-rror of interpolation of positive definite quadratic form by a polyno- 
mial from Pi on all d-dimensional boxes of fixed volume V is 



X d^^2 
A = -Vd. 
4 



(31) 



Proof: Clearly, the minimum of the function with additional assumption (26) is achieved 

when all hi are equal, i.e. hi = h2 = ■ ■ 



i=l 



hd := h. From condition (26) we also have 

d 



, and, hence, A = rnin < = dh'^ = d 



□ 



Now let the quadratic form 



g(x) = J]A 



(32) 



with y4j > for all < i < /c and < for all + 1 < i < c? be given. 



Lemma 7 The L^o - error of interpolation of quadratic form (32) by polynomials Pi on the d- 
dimensional box P of volume V{P) is 



^k-d{d-ky--d I V{P), 



(33) 
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Proof: For the given quadratic form Q{x) = let us consider a linear transformation 

F such that 



i=l 



iQoF)iu) = Y,ul 



i=l 



In other words, 



F(u) 



Ui 



Ud 



Observe that the determinant of the inverse of this transformation is 



(34) 



(35) 



det(F-i) 



(36) 



Let us consider the box F ^{P) which clearly has the volume 

V{F-\P)) = V{P) det(F-^). (37) 
Combining the result of the previous lemma about the error of interpolation on the box F^^{P) 



with (37) and (36), we obtain (33). □ 



Similarly, in the case of positive definite form we obtain the following statement. 

Lemma 8 The error of interpolation of the positive definite quadratic form by polynomials Pi on 
the d-dimensional box P of volume V{P) is equal to 



d 



V{P), 



{3t 



5 Error of interpolation of functions defined on [0, 1]^. 
Estimate from above. 

5.1 Estimate from above for interpolating splines. 

Lemma 9 Let f e Cl{D). Then 

iV^R^(/) 



lim sup 

Af-5-oo 



7M / |i/(/;x)|5fi(x)^rfx 



D 



< 1. 



Proof. 

Let / G C'^{D) be given. For a fixed e G (0, 1) and for every we define 
m^:=min<;m>0: ^ ^* (/, ^) < r , 



(39) 



(40) 
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where uj*{f,S) is a function defined in ([8]). Observe that for m^v defined in such a way it is true 
that m/vT — 7- oo as — 7- oo. In addition, 



— )■ oo, as — )■ oo, 



m 



(41) 



N 



i.e. rriN = o{Nd) as N ^ oo. Indeed, by the definition of niN for all large enough we have 



Nd (rriN — 1) 



m 



■N 



m 



N 



{niN - ly 



> 6 



(rriN - 1' 



m 



N 



2(mN - I] 



oo, as A^ — )■ oo. 



since (^ ""^jv^ ) ~^ ^ ^* 2(ml-i) ) as A^ oo. Hence, (41) is proved. 

Let us divide the unit cube D into cubes with side length equal to ^ and denote the resulting 
cubes by Dj^ , 1 = 1,..., mf^. Next let us take the center point in each cube Df and set 



A. 



In addition, denote by 



N,l 



1 d'f 



xf), i,j = l,...,d, 1 = 1 



2 dxidxj 



, . . . , rrij^ 



HU 



I = 1, 



i=l 



Set the number of elements to be used on Df to be 



n 



N 



Nil-e)\Hi^r)m^. 



I = l,...,m%. 



(42) 



It is essential that nj^ — )■ oo when N oo. This follows from the estimate 



A^(l - £)min{|iJ(/; x)|}2 min{fi(x)} 

1 5 
™(i II fJll 2 II O II 2 
Nil ||oo||ii||oo 



(43) 



together with (41), and the fact that min{r2(x)} > and min{|i7(/; x)|} > 0. 

xS-D xeD 



Now for k ^0, d let us set 



d-k 

1 f d-k^ 2d 



d N I ^ 



l<i<k, 1 = 1, 



1 f d-k\ 2d 



k 



k + 1 < i < d, 1 = 1,... ,m%. 



(44) 
(45) 



12 



For the positive definite form, i.e. A; = or A; = we set 



h 



N _7N 



1 



1 



2 V mfjTif 



I = 1 



, . . . , 



(46) 



(47) 



The intersection of the lattice 

[Lhf,{L + l)h^]''x[Lhf',{L + l)hf']''-\ LeZ, l = l,...,m%, 

with Df^ provides the partition of , I = 1, ...,m^. 

The union of partitions of each Df, I = l,...,m%, provide the partition of D. Denote it by 
□^(D). It can be easily seen that due to definitions (44) and (45), combined with the estimate 
(43), the constructed partition satisfies ([T]) and, hence, is admissible. 

Let us show that the sequence of obtained in such a way partitions {□^(D)}^^^ will be asymp- 
totically optimal. 

By /tv denote the piecewise quadratic function constructed in the following way. On we set 

d 

Jn to be J2 K^A- Then for / > 1 on Df \ Z^f we set 



N,l 2 



i=l 



To estimate the error Hn^J') we observe that 

Rjv(/) < 11/ - S{f, □^)||oo,f^ < 11/ - fN\\oo,n + WfN - s{fN, □^)||oo,Q 

+ MfN, - s{f, D*^)\\oo,n < 211/ - fN\\oo,n + ||/iv - s(/^, □^)||oo,n. 
Let us estimate each term. First of all, by Lemma [T] and the definition of mjv we have 



II/-/ 



Af oo,f7 



< 



1 



2 V2m 



N 



CO* f, 



2m 



N 



I ^ 1 1 oo ^ 2 11^1 



Let us estimate the second term now. Let Rf G \I\*^{D^) be an arbitrary element. By Lemma [7| 
for every x G Rf we have 



|/^(x)-s(/^,n^;x)|fi(xr) < ||/,v-^5(/iv,n^;-)IL^(«-)^(x, 

2 

1 



N\ 



lk4 



d N 



By the definition of nj^ , for all large enough A^, for all I, and for all x G Rf^ , we have 

|/^(x)-s(/^,n^;x)|fi(xf^)< 



< 7k,d 



xf)l^^](xfli 



m%N{l-e)\H{^f)\-^n{^fy. 



fi(xf) 
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Since this estimate does not depend on x, we obtain 



\\fN-s{h,D%)\\^,n< 



7k,d 



Lg|i^(xf)|^^](xf) 



Note that 



mi, 

— ^|iJ(xf)rfi(xf)5^ / |if(/;x)|i(](x)5dx, as N 

^N^^ Jd 



Hence, for all large enough we have 



||/^-s(/^,n; 



Therefore, 



11/ 



N)\\oo,n 



2e 



< 



7k,d 



lk,d 



oo. 



|i7(/;x)|5fi(x)5f/x . 



D 



Ni {N{l-e))i \Jd 



|i/(/;x)|^l](x)^rfx 



Because e > is arbitrary, we obtain the desired estimate from above (39) for RAr(/). 



5.2 Construction of asymptotically optimal sequence of continuous quasi 
interpolating splines. 

In this section, we shall refine the sequence of partitions {nAr}?^=i constructed following the al- 
gorithm in the previous section, to obtain a sequence of partitions and a sequence of continuous 
splines □Ar)}?^=i on the partitions which will interpolate / at all but o{N) as — j- oo points. 

Let parameters of the grid hf and hf^ be as defined in ( |44[ ) and (45) for k ^ 0,d, and in (46) 
for k = 0,d, respectively. Recall that the intersection of the lattice 

[Lhf,{L + l)hf]'x[Lhf',{L + l)W]'-\ LeZ, l = l,...,m%, (48) 

with Df^ provides the partition of Df, I = 1, ...,m^. 

Let us consider two neighboring (i-dimensional regions and with corresponding inter- 
polants and S"^. The partition of consists of parallelepipeds {Dj} and the partition of 
consists of parallelepipeds If the parameters of the grid on and are different, we 

have to subdivide parallelepipeds that have nonempty intersection with the common face to en- 
sure the global continuity of the approximant. Let us assume that the common face lies in the 
{d — l)-dimensional coordinate plane. 

By {Dj^} and {D]i} let us denote parallelepipeds which have nonempty intersection with the 

common face of and D^. In addition, by {Dj^} and {D'^^} we shall denote their {d — 1)- 
dimensional faces contained in fl D^. The set of all possible intersections Dj^ fl constitute 
the partition of fl D^, which is a refinement of each of partitions {Dj^} and 

For each D]~nZ^ we consider D}^ n {(D]~ni^) x M} and n {(D]~ni^) x M}, the set 
of which constitute a refinement of {D}^^} and {D'^i}, respectively. The new vertices of the refined 
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A- 



— ^ 6 — m • 

O Regular vertices on D 1 
• Regular vertices on D2 

A Irregular (new) vertices: interpolate SI (on Dl) or S2 (on D2) 

Common face of Dl and D2 

Additional faces of elements of partition 



Figure 1: Stitching of partitions on neighboring regions in the case d = 2. 



partition we shall call "irregular" to the contrast with vertices of the non refined partitions which 
we shall refer to as "regular" . 

The continuous spline S on two neighboring elements and is constructed now as follows: 
S interpolates at all (regular and irregular) vertices of the partition of D^, and S interpolates S"^ 
at all (regular and irregular) vertices of the partition of D^. Note that automatically the multilinear 
spline S interpolates /, in particular, at all the vertices of the refined partition {D}^^ n -D|i} of the 
common face (1 D"^. 

The illustration of "stitching" of partitions on neighboring regions in 2-dimensional case is shown 
in Figure [1} 

Denote the final partition of D by D*^{D). It can be easily seen that due to (44) and (45), or 
(46), the constructed partition satisfies ([T]). 

Let us compute the number of irregular boxes (denote it by K[^). Recalling (44) and (45) 
together with the fact that the volume of a box from n*j^{Df) is / pf we have that 



K^<c^{ni j^. 

Hence, the total number of irregular boxes is not greater than 



m'^ m' 
C2Y^{nf)-i < C2 

1=1 1=1 



N(l-E)\H(xr)\in{xni 



< CsNdm 



N 



o{N) 



J 



as iV — )■ cx) (since m]\f = o (^^"^ j as — t- oo because of (41)). In particular, this implies that the 
number of boxes in the constructed partition will not exceed N for all N large enough. 
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This verifies the Remark after Theorem 1, because the interpolating spline s{fi\f, D^) constructed 
in the previous section can be discontinuous only along the faces of the irregular boxes. 

It is clear that for two embedded elements of partition the error of multilinear interpolation of 
a quadratic function is greater on the larger element. Therefore, we shall estimate this error on 
elements that do not have intersection with the boundary. Hence, the estima te fr om above for the 
error of quasi interpolation by splines D^) can be obtained as in Section 



5.1 



6 Error of interpolation of functions defined on [0, 1]^. 
Estimate from below. 

Let quantities rriN, nf, Df etc. be as defined in the previous section. 
Lemma 10 Let f G Then 

hminf NiRM) ^ > ^ (49) 

n |i/(/;x)|5n(x)tdx 

Proof. To obtain the estimate from below we shall consider an arbitrary sequence of admissible 
partitions, i.e. box partitions {DArj^^i which satisfies ([T]). 

Note that from (jl| it follows that for an arbitrary element R G Dat, diam(i?) < with some 
constant C > 0. Let us consider the -^^^-neighborhood of the boundary of an arbitrary box Df ^ 
I = 1, ...,m'lf. For an arbitrary e > 0, the volume of the complement of the -^^^-neighborhood of 
the boundary of Dj^ , i.e. the "interior" of D^^ , is 

1 1 / 2Cm^ y^ 1 / 2CmN \ 1-e 

for large enough since mjv = o{N^^'^) as iV — )■ 00. Therefore, for any / = 1, ...,m^, the sum of 
volumes of the boxes which have nonempty intersection with the "interior" of Dj^ (and, due to ([T|, 
lie completely inside of Dj^) is greater than 

Let us show that for any e > and for any large enough there exists index 1^ such that 
the corresponding Dj^ completely contains a "large enough" element Rj^^ G Dat, i.e. element with 
volume greater than —!f-^- Assume to the contrary that there exists eo such that for an arbitrary 



A^o there exists N > Nq such that for all i = 1, ...,m'^ the volume of each box R G Dat having 
nonempty intersection with the "interior" of Df (and, therefore, is completely inside of Df) is less 



than or equal to 



For each i = 1, ...,m^, by uf denote the number of boxes from Dat that are completely inside 

N 



of Df . Note that 



1=1 i=l 

This implies that there exists i* such that uf < nf. Hence, taking into consideration the assumption 



on the volume of each box that is completely inside of Df (or have nonempty intersection with the 
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"interior" of -DiY) to be less than or equal to % , the total volume of all such boxes from Dat is 
less than or equal to 



m 



N 



which contradicts to (50). 



For each N and corresponding Z^r, set 

d 

/iv,i^(x) := ^Af/xl 



i=l 



Observe that 

By Lemma [7] we have for some e > 

1 



\\fN,lr, - s(/iV,ijv, n^)llLoo,n(i?r ) ^ ~ ^)^M 



rf IV 



By the definition of nf we have that for all large enough 



lk,d 



> 



IkA 



Ik d 
Nd 



i=i 

m^Af{l-£)|i?(xf')pn(xf)2 



d V ,^(x, 



v 

^|/7(/;x)|i^](x)irfx 



Hence, for all large enough we obtain 



>7m 



|i/(/;x)|^fi(x)5dx 



On the other hand 



11/ - /w,iivllLoo,n(fi,^) ^ 11^"'^]^ 



due to the choice of 111^. Hence, we obtain that for all large enough N 



||/-s(/,n;v)||oo,n> (1-C4e)^ 



D 



with some positive constant C4. Therefore, 



liminf 11/ - s(/, □iv)||oo,n > _2 



Af-s>oo 



7M 



|i/(/;x)|2(](x)2cix . 



D 



□ 



The estimate from below for quasi interpolating sphnes can be obtained analogously. 
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